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Effect of Internal Radiant Heat Transfer on Temperature 
Distribution, Thermal Stress, and Deflection in Box Beams l 

Stanley Goodman, Stanton B. Russell, and Charles E. Noble 

Thirteen box beams were analyzed to determine the effect of internal radiant heat 
transfer on temperature distribution, thermal stress, and deflection. A rough experimental 
check of temperature distribution and beam deflection was made for one case. For maxi- 
mum beam temperatures above 700° to 900° F, change in omittance of the interior surfaces 
of the beam had an appreciable effect on the cover (but not the web) temperatures and to 
an even greater extent, on the beam deflection. At maximum beam temperatures of 1,200° F, 
an increase in interior surface emittance appreciably relieved the maximum thermal stress. 



1. Introduction 

The subjection of aircraft and other structures to 
increasingly high temperatures gives rise to a need 
for predicting the transient temperature distribution 
in the structures, and the resulting thermal stresses 
and deformations. This knowledge will make pos- 
sible the rational design of the structure, and the 
devising of methods of alleviating thermal stresses 
and deflections. 

As the temperature of a structure increases, radia- 
tion becomes more important as a mode of heat 
transfer [l]. 2 The influence of radiant heat transfer 
on the temperature distribution of the structure is 
affected by structural geometry, heating rate, and 
thermal properties of the beam material. The 
purpose of this investigation was to determine the 
effect of changing the emittance of the interior sur- 
faces of a box beam on temperature distribution, 
maximum thermal stress, and beam deflection in 
heated box beams. 

2. Structures and Material 

Thirteen analyses were made of modifications of 
three basic types of box beams. The beam cross 
sections are shown in the figure at the top of table 1 ; 
wall thicknesses are indicated by the ratios listed in 
columns 2 and 3 of the table. One type (cases 1, 2, 
and 3) is a thick-walled beam, similar to a beam used 
in the experimental work. The other types are typi- 
cal of beams used in aircraft structures; one type is 
relatively thick-walled (cases 4 to 8) and the other 
relatively thin-walled (cases 9 to 13). 

Thermal and elastic properties were taken as 
those of the type 302 stainless-steel test beam. Ther- 
mal conductivity and specific heat were approxi- 
mated by the linear relationships : 



&=7.08 + 0.0043r 

c = 0.106 + 0.0000257r, 



(1) 

(2) 



1 This work was conducted at the National Bureau of Standards under the 
sponsorship and with the financial assistance of the Office of Naval Research. 

2 Figures in brackets indicate the literature references at the end cf this paper. 



where k is thermal conductivity (Btu/hr ft °F); 
c is specific heat (Btu/lb ° F); and 7 7 is absolute tem- 
perature (°R). Hemispherical total emittance was 
taken as 0, 0.35 (emittance of the test beam, cases 
1 to 3) and 1, as indicated in columns 4 and 5 of 
table 1 . 

Young's modulus and the coefficient of linear 
thermal expansion were represented by 



E= 3 1.0 X10 6 - 0.0064 7 T X10 6 
a=8.29X10- 6 + 0.00137TX10- 



(3) 

(4) 



where E is Young's modulus (lb/in. 2 ), a is coefficient 
of linear thermal expansion (in. /in. ° F), and T is 
absolute temperature (° R). 

3. Boundary Conditions 

The beams were assumed to be subjected to a 
uniform heat input along one cover in a rarefied 
atmosphere. The heating rates are listed in columns 
6 and 7 of table 1. The initial beam temperature 
was assumed to be 78° F. For the test, beam 
(cases 1, 2, 3), the boundary condition was taken 
as the heated cover temperature (figs. 2 and 7) 
observed in an experiment. For all other cases, the 
boundary condition was taken as the assumed heat 
input to the cover. All elements of the heated 
cover were assumed to remain constant at 1,200° F 
after reaching that temperature. 

4. Method of Analysis 

4. 1 . Calculation of Temperature Distribution 

Analysis was made by a numerical method similar 
to that described by Dusinberre [2]. The beam 
cross section was divided into analysis elements, 
and a heat-balance equation was set up for each 
analysis element and solved for the element tem- 
perature after a short time interval, A??. The 
process was repeated for successive short-time 
intervals using the new element temperatures as 
the starting point of the next computation. It 
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Table 1. Computed cover temperatures , maximum thermal stress and maximum deflection in box beams heated along one cover 



\ = 4.85fo -* 



4.85 in. 



1.76 in. 



Case 


ri 


T2 


Surface emittance 


Heating rate 


Elapsed 
time 


Temperature b at 


Maximum stress c 


Maximum 
deflection c d 


Interior 


Exterior 




A 


B 


Tensile 


Compressive 


1 
2 
3 


0. 0727 
.0727 
.0727 


0. 0264 
.0264 
.0264 




0.35 

1 


0.35 
.35 
.35 


BtuJhrfP 


°Fjsec a 
3.8 
3.8 
3.8 


sec 
378 
378 
378 


o F 

958 
958 
958 


O p 

84 
147 
377 


psi 
35. 3X10 3 
35.8 
34.8 


psi 
29.3X10 3 
29.7 
30.5 


in. /in. 2 
0. 46X10-3 
.42 
.30 


4 
5 
6 


.02 
.02 
.02 


.03 
.03 
.03 




0.35 

1 


.35 
.35 
.35 


32400 
32400 
32400 


50 
50 
50 


26.6 
26.6 
26.6 


1202 
1173 
1095 


78 
101 
188 


121.7 

118.0 
108.5 


127.9 
124.4 
116.1 


.49 
.47 
.41 


7 
8 


.02 
.02 


.03 
.03 


0.35 
.35 



1 


32400 
32400 


50 
50 


26.6 
26.6 


1200 
1104 


102 

98 


121.5 
111.0 


127.1 
116.7 


.48 
.44 


9 
10 
11 


.01 
.01 
.01 


.01 
.01 
.01 




0.35 

1 


0.35 
.35 
.35 


10800 
10800 
10800 


30 
30 
30 


45 
45 
45 


1204 
1135 
959 


78 
154 
371 


102.9 
92.1 
71.5 


99.0 
91.0 
74.9 


.58 
.51 
.34 


12 
13 


.01 
.01 


.01 
.01 



1 


.35 
.35 


32400 
32400 


91 
91 


13.7 
13.7 


1200 
1125 


78 
203 


127.6 
114.5 


101.1 
95.1 


.58 
.49 



a Approximate temperature rise rate at point A up to about 800° F. 

b Initial temperature 78° F. 

c Neglecting yielding, creep, and buckling. 

d Deflection per unit beam length squared. 



was assumed that surface emittance remained 
constant and that beam surfaces reflected diffusely 
and were gray, i.e., total emittance and absorptance 
are equal at all beam temperatures. Gaseous heat 
transfer was neglected. 

One-half of the cross section of the symmetrical 
beam was divided into 12 analysis elements. Since 



the external heat transfer was symmetrical with 
respect to the left and right halves of the beam 
cross section, the right half of the beam was repre- 
sented by a perfect but diffuse reflector. The 
reflector was divided into four regions. The element 
configuration for cases 1, 2, and 3, shown in figure 1, 
is typical of all cases. 



10 II 



2.02' 



[^ 



12 



^ 



14 

-I- 



15 



16 



«*- 0.128 



h-0.5ll 

— 2.05 

2.472" 



'4«— 1.022" J 



W- 



JJ 



Figure 1. Test beam cross section (cases 1. 2 1 3). 

Material: Stainless steel, 18-8 (302); length: 30 in.; analysis elements designated by dotted lines and 
small numerals. 
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The heat balance equation is as follows : 



T' = T 



Ad 



;{ 



\J€m—l,m)\-H- m )\l m — i 1 m ) 



V mP^m 

+<rA^[F m , 1 (Tl-T,i)+F m . 2 (TI-Ti)+ . . . 
+^, 12 (? 7 1 4 2-^)]+T^\, 



where 



}• 



T' m is the temperature (°R) at the center of element 
m after time interval Ad 

V m is the volume (ft 3 ) of element m 

T m is the temperature (° R) at the center of 
element m at the beginning of the time interval 
Ad 

p is the density of the material (lb/ft 3 ) 

C m is the specific heat of element m (Btu/lb ° F), 
represented by a linear function of T m . 

k m , m+1 , k m - lf m> are the thermal conductivities of 
the beam metal at the mean temperatures of the 
elements m and m+1, and m and m—1, 
respectively (Btu/hr ft ° F) 

H m , J m are ratios of element contact areas to 
distances between element centers for element 
m and adjacent elements m—1, m + 1, respec- 
tively (ft) 

t m is hemispherical total emittance of element 
surface m 

ex is the Stefan-Boltzmann constant (0.1713X10 -8 ) 
in Btu/ft 2 hr ° R 4 

A m is exterior surface area (ft 2 ) of element m 

T A is ambient temperature (° R) 

A' m is interior surface area (ft 2 ) of element m 

F m<n is an over-all radiant heat interchange factor 
for net radiant heat exchange between a gray 
surface of element m and a gray surface of 
element n. It includes the effect of direct and 
all reflected radiation. 

W m is the external heat input (Btu/hr). 

The first two terms inside the brackets on the 
right-hand side of eq (5) represent conductive heat 
transfer between element m and adjacent elements; 
the third term represents radiant heat transfer with 
the outside environment, the fourth term represents 
internal radiant heat exchange between element m 
and the rest of the beam interior; the fifth term 
represents the assumed heat input boundary condi- 
tion. For eases 1, 2, and 3 in which the boundary 
condition was taken as the observed heated cover 
temperature history of the test beam, the temper- 
ature of each heated cover element was expressed by 
four linear functions of time covering successive 
time intervals. 

To compute the graybody radiant heat exchange 
factors, F, blackbody radiant heat exchange factors 
were computed for all two-element combinations of 
the 16 internal surfaces of the analysis elements and 
(fictitious) reflector using the methods of reference 
[3], and radiation geometry. Values of the 256 F's 



were then obtained from the matrix relationships of 
reference [4]. The computations were performed on 
SEAC using an existing code for inverting the 
mat rices. 

The temperature-distribution computations were 
computed on SEAC. Temperatures at the mid- 
points of the 12 analysis elements of the beam were 
printed out at regular time intervals which were 
whole number multiples of Ad. 

In order to minimize computing machine running 
time, it was desirable to use the largest value of the 
time interval, A0, consistent with adequate accuracy 
of solution. The time interval must be sufficiently 
short that during the interval (1) The initial element 
temperatures can be used with negligible error and 
(2) for conductive heat transfer; the effect of any 
nonadjacent element on a given element is negligible. 
Satisfactory values of Ad were obtained by trial. A 
portion of the temperature distribution history was 
computed several times using successively larger 
values of Ad. The largest value of Ad which gave a 
temperature distribution history negligibly different 
from that obtained with the smallest value of A0 was 
used to compute the complete temperature distribu- 
tion history. Values of Ad found satisfactory by this 
method ranged from 0.72 sec for rapidly heated 
beams to 1.8 sec for slowly heated beams. 

Some error was generated in the beam temperature 
history computations by representing one-half of the 
beam cross section by a perfect, diffuse rather than 
specular, reflector. As a check on the magnitude of 
this error, radiant heat transfer rates were computed 
for a simple symmetrical case using first the entire 
beam cross section in the computations, and then an 
equivalent beam consisting of one-half of the beam 
cross section and a diffuse reflector. An infinitely 
long box beam of rectangular cross section divided 
into six analysis elements was used. The beam was 
2-in. deep, 8-in. wide, and emissivity was 0.35. The 
vertical (2-in.) walls of the beam were at temper- 
ature Ti and the horizontal walls at absolute zero. 
The rate of radiant heat transfer to the left half of 
the upper horizontal wall was 

#=O.13540-jTi for original beam 

^— 0.1 389(r T\ for equivalent beam 

with diffuse reflector, (6) 

where q is rate of radiant heat transfer (Btu/hr ft 2 ); 
a is the Stefan-Boltzmann constant (0.1713 X10~ 8 
Btu/(ft 2 hr °R 4 )); and Ti is vertical wall tempera- 
ture (° R). It was concluded that the error due to 
use of a diffuse rather than a specular reflector was 
small. 



4.2. Computation of Thermal Stress and Deflection 

Thermal stresses and deflections in the beam were 
computed by the methods described in reference [5]. 
Integrations over the beam cross section were 
performed numerically using finite elements identical 
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with those used in the temperature distribution 
computations (fig. 1). Possible effects of yielding, 
creep, and buckling were neglected. 

5. Test of Beam 

The test specimen, corresponding to case 2, table 
1, was a box beam constructed of 0.128-in. -thick 
type 302 stainless-steel sheet (fig. 1). Two pieces 
of sheet were each bent into identical channel shapes 
and butt welded together lengthwise to form the 
2 in. by 5 in. by 30 in. beam. 

The beam was heated in a vacuum chamber whose 
pressure w^as maintained at 4 mm of mercury, a 
pressure low enough to eliminate appreciable con- 
vective heat transfer. Only one cover of the beam 
was subjected to heating. It was heated with 12 
quartz-tube, tungsten-filament heating elements uni- 
formly spaced in a reflector. Total output of the 
heaters was 5.7 Btu/sec. The reflector, a rectangular 
silver-plated stainless-steel channel enclosing the 
heating elements, confronted the heated cover of the 
beam. The reflectance of the reflector surface was 
about 0.97. 

Temperatures were measured at the center cross 
section on the heated cover of the beam, at the 
center cross section on the unheated cover, at two 
points on the heated cover 1 in. from one end, and at 
one point on the heated cover 1 in. from the other 
end. The temperatures were measured with iron- 
constantan thermocouples whose outputs were indi- 
cated by galvanometers. The thermocouples were 
calibrated before and after the test. The temper- 
ature rise of the heated cover was approximately 12 
percent less at the beam ends than at the center 
cross section. The experimentally determined heated 
cover temperatures, which were used as the boundary 
condition in cases 1,2, and 3, are shown in figure 2 
and, together with observed temperatures for two 
points on the opposite cover, in figure 7. 

The normal total emittance of the type 302 
stainless steel was determined by comparing its 
rate of radiation at a given temperature with that of 
a silicon carbide bar at the same temperature. The 
radiation rates were determined by focusing the 
images of equal areas of the silicon carbide bar and of 
the stainless steel successively on a thermopile, using 
a fluorite lens. The emittance of the stainless steel 
obtained was constant in the temperature range 400° 
to 1,000° F and equal to 0.35. 

Beam deflection at the center cross section was 
measured by means of SR-4 type AB-5 electrical 
strain gages mounted back-to-back on a shielded 
cantilever beam which was deflected by displace- 
ment of the center of the beam relative to its ends. 
Contact of the center of the beam with the cantilever 
was made by a Vycor rod. Calibration was ac- 
complished by deflecting the cantilever by known 
amounts at the point of contact with a micrometer 
screw. 

The maximum deflection observed was 0.415 in. 
after 240 sec of heating. After 378 sec of heating 
the center deflection was 0.37 in. 



6. Results 

6.1. Beam-Temperature Distribution 

Computed temperature histories at point A, in 
the center of the heated cover and at point B, in the 
center of the opposite cover (see sketch on table 1) 
are shown in figures 2 to 6 for the 13 cases considered. 
Figure 7 shows the computed temperature distribu- 
tion in a thick-walled beam (cases 1, 2, 3) of the 
same dimensions as the test beam after 148, 189, and 
297 sec of heating. Figure 8 shows the temperature 
distribution in a thin-walled beam after 45 sec of 
slow heating, and after 13.7 sec of fast heating. 
Temperatures at A and B are listed in columns 9 
and 10 of table 1. Comparison was made of geo- 
metrically similar beams heated at the same rate 
for the same length of time but having the following 
different values of interior or exterior emittances: 0, 
0.35, and 1.0. 

It was found that, in all cases considered, a change 
in emittance of the surfaces had an appreciable effect 
on the cover temperatures and little on the web 
temperatures when the maximum heated cover 
temperature was above 700° to 900° F. For beams 
with the heated cover at a temperature of about 
1,200° F, changing the internal emittance from to 
1.0 reduced the temperature difference, d, between 
points A and B by 18 percent for the thin-walled 
beam with the cover heated at 91° F/sec, and by 48 
percent for the same beam with one cover heated at 
30° F/sec. For the same cover-heating rate the 
effect of change in interior emittance was less for 
the thicker-walled beam. 
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Figuee 2. Beam temperature history. 
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Figure 3. Beam, temperature history for heating rate of 3%, 400 
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Figure 4. Beam temperature history for heating rate of 32,400 
Btujhr ft 2 . 
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Figure 5. Beam temperature history for heating rate of 10,800 
Btujhr ft 2 ." 
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Figure 6. Beam temperature history for heating rate of 32,4.00 
Btu/hr ft 2 . 
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A change in exterior surface emittance has less 
effect on d than a change in interior surface emit- 
tance. Comparison of cases 4 and 6 with corres- 
ponding cases 7 and 8 indicates that, for the particu- 
lar conditions specified in table 1, increasing the 



interior surface emittance from to 1.0 causes a 
decrease in d of 19 percent; increasing the exterior 
surface emittance from to 1.0 causes a decrease in 
d of only 8 percent. 
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Figure 7. Test beam temperature distribution. 
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6.2. Thermal Stress and Deflection 

In order to evaluate the effect of a change in 
emittance on the thermal stress distribution and 
deflection, thermal stress distributions and deflec- 
tions were computed for the 13 cases under the 
conditions listed in table 1. Maximum tensile 
stress, maximum compressive stress, and maximum 
deflection are listed in columns 11, 12, and 13, 
respectively, in the table. To generalize the results 
as regards beam length the deflection is given as the 
deflection per unit beam length squared; the maxi- 
mum deflection is proportional to the beam length 
squared, the thermal stresses are independent of the 
beam length. The possible effects of yielding, creep, 
and buckling were neglected. Values of the thermal 
stress and deflection are therefore valid only for 
comparison purposes. 

For a thick-walled, slowly heated beam with 
heated cover at 958° F (cases 1, 2, 3) a change in the 
interior wall emittance has little effect on the magni- 
tude of the maximum thermal stress. An increase 
in the interior wall emittance from to 1.0 however, 
results in a decrease of 34 percent in the maximum 
beam deflection. For the beams with the same 
geometry and under the same heating conditions in 
which the final heated cover temperature is about 
1,200° F, a change of the interior wall emittance 
from to 1.0 reduces the absolute maximum stress 
by from 9 percent (cases 4 and 6) to 27 percent (cases 
9 and 11), and deflection by from 15 percent (cases 
12 and 13) to 40 percent (cases 9 and 11). In the 
latter beam, the reduction in maximum tensile stress 
was 31 percent. For a change of the exterior wall 
emittance from to 1.0 (cases 7 and 8), the reduc- 
tion in absolute maximum stress was 8 percent, and 
the reduction in maximum deflection was 7 percent. 

For the beams analyzed, changing the interior 
beam emittance from 0.35 to 1 causes a change in 
the temperature difference, d, between A and B 
which is about 75 percent of the change that would 
occur for a change in emittance from to 1. Cor- 
responding ratios for maximum thermal stress and 
deflection are about 66 and 74 percent, respectively. 



6.3. Comparison of Experimental and 
Theoretical Results 

A partial check of the temperature distribution 
and deflection was made for case 2 using a beam 30 
in. long. Computed and observed temperatures for 
two points on the unheated cover after 148, 189, 
and 297 sees of heating are shown in figure 7. Agree- 
ment is fair. Some of the discrepancy may be 
attributed to gas heat transfer within the beam. 
For the conditions listed under case 2 in table 1, 
observed beam deflection was 0.37 in. and beam 
deflection computed from the theoretical temper- 
ature distribution was 0.38 in. i.e., 0.00042 in. /(in. 
length) 2 . 

7. Conclusions 

A change in the interior surface emittance of a 
box beam heated along one cover has an appreciable 
effect on the temperature distribution, and con- 
sequently on the maximum thermal stress and the 
deflection. For a variety of beam geometries and 
heating conditions, a change of the interior surface 
emittance from 0.35 to 1 for a stainless steel beam 
reduces considerably the temperature difference 
from cover to cover in the temperature range above 
700° to 900° F and the maximum thermal stress 
and the deflection in the temperature range above 
900° to 1,200° F. The effect is particularly marked 
for a thin-walled slowly-heated beam. 
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Important Notice 

Beginning in July 1959, the Journal of Research will be published in four 
separate sections which may be subscribed for individually. This change is being 
made to provide a publication whose subject matter concentrates on particular 
fields of interest, thus meeting the specialized needs of individual scientists, 
engineers, and mathematicians. 

Section A. Physics and Chemistry, to be issued bimonthly, will present 
papers of interest primarily to scientists working in these fields. This section 
will cover a broad range of physical and chemical research, with major emphasis on 
standards of physical measurement, fundamental constants, and properties of 
matter. 

Section B. Mathematics and Mathematical Physics, to be issued 
quarterly, will present results in pure and applied mathematics, including mathe- 
matical statistics, theory of experiment design, numerical analysis, and short 
numerical tables; theoretical physics, chemistry, and engineering, where the 
emphasis is on the mathematical content or methodology; and logical design, 
programing, and applications of electronic computers and computer systems. 

Section C. Engineering and Instrumentation, to be issued quarterly, 
will report results of interest chiefly to the engineer or applied scientist. This 
section will include many of the new developments in instrumentation resulting 
from the Bureau's work in physical measurement, data processing, and the develop- 
ment of test methods. It will also cover some of the work in acoustics, applied 
mechanics, building research, and cryogenic engineering. 

Section D. Radio Propagation, to be issued bimonthly, will report research 
in radio propagation, communications, and upper atmospheric physics. Topics 
to be covered include propagation in ionized media, scattering by turbulence, the 
effect of irregular terrain on propagation, diffraction and scattering by solid 
obstacles, propagation through time-varying media, and antennas. 



All present subscribers to the Journal will receive an announcement from the 
Superintendent of Documents when advance orders for the new sections can be 
accepted. 
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